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^The  basic  objective  of  this  experimental  study  Is  to  optimize  the  process 
of  Injecting  drag-reducing  additives  into  a  water  flow.  The  Initial  phase 
consisted  of  a  factorial  design  where  the  concentration  and  flowrate  of  the 
Injected  additive  as  well  as  the  angle  and  width  of  the  Injection  slot  were 
varied  independently.  The  experiments  were  conducted  In  a  rectangular  cross 
section  channel  that  has  an  aspect  ratio  of  ten  to  one.  Slots  were  located  in 
both  of  the  larger  walls  at  a  streanmlse  location  where  the  channel  flow  of 
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^>water  was  fully  developed.  Drag  reduction  was  deduced  from  wall  pressure 
measurements  3nd  wail-layer  concentration  of  the  additive  was  deduced  from 
colorimeter  measurements.  Additive  concentration  varied  from  100  to  400  ppm. 
Injection  flowrates  ranged  from  200  to  400  ml/min.  The  slot  angle  was  either 
15  or  25  degrees  while  the  slot  width  was  either  1.27  ran  or  2.54  mm.  The 
quantity  which  was  optimized  (the  merit  function)  was  the  integral  of  the  drag 
reduction  with  respect  to  streamwise  distance  beginning  at  the  initial  location 
of  positive  drag  reduction  and  ending  at  the  location  where  the  additive  was 
uniformly  mixed  with  the  water. \ 


Statistical  analysis  of  the  results  showed  that  for  the  range  of  variable 
studied,  the  merit  function  depends  primarily  upon  the  linear  terms  for 
additive  concentration  and  injection  flow  rate.  The  product  of  additive 
concentration  and  injection  flowrate  as  well  as  the  product  of  slot  angle  and 
slot  width  were  significant  but  less  important  factors  in  fitting  the  experi¬ 
mental  response.  Tne  analysis  indicated  that  the  largest  Increase  In  the 
merit  function  would  be  achieved  by  increasing  the  injection  flow  rate  and 
additive  concentration. 

The  results  showed  chat  additive  concentrations  less  than  1  ppm  yield  drag 
reduction  on  the  order  of  20%.  These  results  were  achieved  well  downstream  of 
the  Injector  where  it  is  hypothesized  that  the  additive  molecules  were  in  an 
excellent  conformation  for  reducing  drag.^_ 
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INTRODUCTION 


It  1*  veil  known  that  the  presence  of  dilute  solutions  of 
fcifli  aoiicular  weight  polyaers  in  tho  near-well  region  of  bound 
turbulent  shear  flows  can  reduce  the  wall  shear  stress.  As  shewn 
by  recent  studies  at  Purdue  (1*2)  and  by  HcCoab  and  Rabie  <3)«, 
the  additives  have  a  direct  effect  upon  the  flow  structures  in 
the  buffer  lager ,  10  <  y*  <  100.  The  purpose  of  this  study  was 
to  develop  eiperia^ntally*  optlaua  aetheds  for  the  injection  of 
polyaer  additives  into  a  turbulent  water  flow.  Tho  specific 
objectives  were  to  aaslaizo  tho  reduction  of  viscous  drag  down** 
streaei  of  the  injection  slot  and  to  build  a  data  baso  for  the 
diffusion  characteristics  of  tho  injected  additives. 

These  Oiperlaents  were  conducted  in  a  2.  9  ca  by  29  ca  rec¬ 
tangular  cross  section  channel.  The  additive  solutions  were 
injected  into  the  fully  developed  turbulent  water  flow  through 
slots  in  the  29  ca  walls.  Drag  reduction  was  deduced  froa  pres¬ 
sure  drop  aeasureaents  aade  upstreaa  and  downstreaa  of  tho  injec¬ 
tion  slots  with  and  without  additive  injection.  Heaswreaents  of 
the  additive  concentration  in  the  near-wall  region  wore  obtained 
by  withdrawing  saaples  of  fluid  froa  this  region  during  injec¬ 
tion. 

The  eptlaiiation  progrea  wee  begun  by  testing  all  c sab 1 na¬ 
tions  of  two  levels  of  each  independent  variable  governing  the 
injection  and  diffusion  processes.  These  were  the  angle  and 
width  of  the  injection  slot  and  the  concentration  and  flowrate  of 


the  additive  solution.  Tho  results  of  these  experiments  were 
analysed  statistically  using  the  techniques  of  response  surface 
Methodology  to  determine  the  most  efficient  way  to  approach  the 
optimum  combination  of  those  four  variables. 

There  have  been  several  other  experimental  programs  where 
the  injection  process  has  been  studied.  For  example  Hsus  and 
Wilhelm  (4)  examined  the  effect  of  additive  injection  on  drag  in 
a  pipe  with  five  circumferential  slots  located  six  inches  apart. 
While  they  varied  injection  flowrates  and  determined  an  optimum 
injection  concentration  for  each  injection  pattern*  they  did  not 
vary  the  design  of  the  slots. 

Walters  and  Wells  (9)  discussed  the  desirability  of  reducing 
an  anomalous  drag  Increase  detected  in  their  studies  of  injection 
through  a  porous  wail.  They  varied  injection  flowrate  and  injec¬ 
tion  concentration  for  different  lengths  of  porous  wall  injec¬ 
tion.  However*  similar  to  Haus  l>  William  there  was  not  an 
independent  variation  in  the  porosity  of  the  injector. 

Hore  recently  an  anomalous  drag  increase  near  the  injector 
was  reported  and  discussed  by  Fruman  and  Oelivel  (6*7).  In  addi¬ 
tion  to  decreases  in  mixing  as  concentration  increases  they 
reported  swelling  of  the  injected  material  as  it  leaves  the  slot. 

Th w  studies  of  Wu  and  Tulin  (8)  and  Wu  (9)  included  sys¬ 
tematic  variations  in  injection  flowrate*  concentration  of 
injected  fluid  and  injector  design.  However*  in  those  two  stu¬ 
dies  the  total  drag  on  a  flat  plate  of  fixed  streesnxise  length 
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was  measured.  In  this  study  the  static  imU  pressure  downstream 
of  the  injection  was  Measured.  Froa  th i <  streaawise  pressure 
distribution  one  directly  measures  any  adverse  dray  increase  and 
determines  the  streemwise  length  over  which  the  injection  is 
effective.  The  latter  is  crucial  for  determining  the  need  and 
location  for  a  second  slot. 

The  present  experiments  were  conducted  in  the  fully 
developed  region  of  a  turbulent  channel  flow.  Consequently  they 
differed  significantly  from  the  studies  at  Colorado  State  (10*11) 
where  the  drag-reducing  additive  was  injected  into  the  strong 
favorable  pressure  gradient  in  the  entrance  region  of  a  plpo. 

There  have  also  been  several  studies  where  the  diffusion  of 
a  'drag-reducing  additive  from  the  wall  region  has  been  studied 
( 12* 13.  14.  13).  While  these  studies  have  provided  valuable 
insights*  including  the  conclusion  that  large  injection  concen¬ 
trations  can  yield  unnecessarily  large  concentrations  in  the  wall 
region  (14).  they  have  not  yet  lead  to  methods  for  achieving  an' 
optimum  design  for  the  injection  process. 
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EXPERIMENTAL  APPARATUS  AND  PROCEDURES 

Aanratm 

The  experlwent*  vara  perforated  in  the  recirculating  flow 
loop  shown  in  Figure  1.  The  flow  loop  incorporates  a  c oath inet ion 
of  e  perforated  plate*  a  screen-sponge-screen  section*  saooth 
contractions  and  flow  straighteners  to  ensure  a  swoo«h  flow  at 
the  inlet  of  the  two  diatensional  channel.  The  flowrate  in  the 
channel  is  aionitored  using  an  orifice  water  connected  to  an 
inverted  U-tube  ownoaeter. 

The  test  section  of  the  channel  has  a  rectangular  cross-* 
section*  with  internal  diwensions  of  2.9  cw  by  29  'em  resulting  in 
an  aspect  ratio  of  ten  to  one.  The  injection  slots  are  located 
wore  than  one  hundred  channel  widths  downstreaa  of  the  inlet  and 
approxiaMtely  eighty  channel  widths  upstreaw  of  the  exit.  Conse¬ 
quently*  the  flow  at  the  injection  point  is  typical  of  fully 
developed*  two-diwenslonal  channel  flow. 

The  additive  solutions  were  injected  through  reawyahle  slots 
located  in  both  of  the  29  cw  walls  of  the  tost  section.  Pressure 
drop  weaeurewonts  woro  obtained  frow  thirteen  proeeuro  tape 
located  along  the  centerline  of  one  of  the  29  cw  walls.  The 
relative  locations  of  the  slots  and  pressure  tape  are  shown  in 
Figure  2.  The  injection  elate*  alto  shown  in  Figure  2*  are  in 
braes  inserts  that  span  the  channel  in  the  29  cw  walls.  The 
slots  have  a  length*  in  the  epanwiee  direction*  of  22. 44  cw  and 
their  geowetry  is  defined  by  the  slot  width  (d)  weasured  in  the 


Schematic  of  flow  loop. 


-  7  - 


atreaawiae  dirtciion  and  the  angle  of  inclination  of  th#  outlet 
to  the  flow  direction  in  tha  channel  <a> .  Th#  slots  are  aechined 
with  an  included  angla  of  five  degreea  and  tha  angle  of  inclina¬ 
tion  ia  Matured  to  tha  bleactor  of  thia  includaa  angla.  For 
this  study,  slot  width*  of  2.94  mm  and  t.  27  mm  were  uaod  along 
with  alot  anglaa  of  19  and  29  dagraaa. 


Praaaura  drop  aaaaureaent*  wara  aada  with  two  Olluont 
aicroaeter  aanonatara.  Tha  aanoaeter  fluid  utad  waa  carbon 
tetrachloride  raaultlng  in  a  aaaaureaent  eenaitivlty  of  0.019  mm 
of  watar. 


Yha  additlvaa  ware  diluta  aolutiona  of  SEP ARAN  AP-273.  a 
polyacry 1 laid#  aanufactured  by  Dow  Cheaical  Corp.  Solutiona 
with  concentratlona  of  100  ppw  and  400  ppa  by  weight  wara  gravity 
fad  to  tha  injector  alota  fro*  an  ovarhaad  raaarvoir.  Tha  injec¬ 
tion  flowrate  waa  regulated  uaitig  a  watering  valve  And  a  Oilaont 
rotaaatar  flow  aatar  for  each  alot.  The  flow  aatjire  wara  cali- 
eratad  for  each  concentration  of  tha  additive. 

Concentration  aaaauraaanta  wara  aada  uaing  a  fauach  and  !Loah 


8pactronlc  20  8pectropfcoto*eter  to  aaaaura  tha  dyo 


concentration 


in  aaaploa  drawn  frtw  the  channel  while  Injecting  dyad  additive 


kolutiona.  Uaing  an  initial  dye  concentration 
Ploureaeein  dlaodiua  aalt  raaultad  in  an  uncertainty 
odd*  of  four  percent  for  aaaaurad  dye  concentratlona 
one  percent  of  tha  initial  dyo  concentration.  Tha 


of  2  g/1  of 
at  20  to  1 
greater  than 
uncertainty 


lncraaaad  to  twenty  percent  for  aaaaurad  dye  ebneentratione 
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between  one-tenth  of  one  percent  end  one  percent  of  the  initial 
dye  concentration. 

Procedure* 

Prior  to.  an  experiment,  filtered  softened  tap  water  for  the 
flow  loop  was  deaerated  by  heating  to  approximately  50°C  in  a 
separate  holding  tank.  This  water  war  then  cooled  to  room  tem- 
perature  before  being  introduced  into  She  flow  loop. 

Preparation  of  the  additive  solutions  consisted  of  two 
steps.  The  solutions  were  initially  mixed  to  a  concentration  of 
2670  ppm  and  allowed  to  hydrate  for  twelve  to  twenty-four  hours. 
These  solutions  were  then  diluted  to  their  final  concentration, 
either  100  ppm  or  400  ppm.  and  allowed  to  hydrate  for  another 
twelve  to  twenty-four  hours  before  boing  used.  Filtered  tap 
water  was  used  in  the  preparation  of  all  solutions. 

Prior  to  each  experiment  the  drag  reducing  capability  of 
each  batch  of  additive  was  established  using  a  horizontal  1.409 
cm  X. 0.  tube.  The  additive  was  gravity  fed  to  the  tube  with  the 
flowrate  controlled  by  a  valve  at  the  tube  outlet.  The  flowrate 
through  the  tube  was  measured  along  with  the  pressure  drop  across 
two  taps,  one  located  more  than  30  tube  diameters  downstream  of 
the  inlet  and  the  other  more  than  19  disasters  upetream  of  tha 
outlat.  Tha  distorts  between  these  tape  ie  two  meters.  The 
viscosities  of  the  additives  were  aeesured  with  a  Wells* 
Brookfield  LVT-SCP  l."t65°  cone  and  plate  micro  viscometer  at 

sec*1 


shear  rates  of  115  and  230 
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Drag  reduction  was  calculated  frow  pressure  drop  ataiurt- 
oents  taken  with  and  without  additive  injection.  The  sequence  of 
•tep*  for  the  preeeure  drop  neasuroMents  were  ae  follows.  First 
with  the  poops  turned  off  and  water  in  the  channel*  a  taro  read¬ 
ing  for  the  aaneaeter  was  taken.  The  loop  was  then  started  and 
the  pressure  drop  across  a  pair  of  taps  was  Measured.  Once  the 
"without-injection"  pressure  drop  was  Measured*  the  injectors 
were  turned  on  and  the  pressure  drop  with  additive  injection  was 
Measured.  The  injection  was  then  ter  'inated  and  the  "wlthout- 
injectlon"  pressure  drop  was  rsMet sured.  Finally  the  flow  was 
stopped  and  the  zero  was  reweasured.  If  the  second  "wlthout- 
injectlon"  pressure  drop  and  zero  oeasureoents  were  In  agreeoent 
with  the  first  ones*  the  Manoweter  was  connected  to  the  next  pair 
of  pressure  taps  in  preparation  for  *  e  next  sequence  of  oeasure- 
aents.  Xn  the  event  that  the  repeated  readings  did  net  reproduce 
the  initial  aeasureMents*  the  sequence  was  repeated  for  that  pair 
of  taps. 

81nce  these  eiperlaents  were,  conducted  in  a  recirculating 
flow  loop*  there  is  the  potential  for  an  aceuswlation  of  additive 
in  the  water*  resulting  in  drag  reduction  without  injection. 
Therefore*  the  condition  of  the  water  in  the  flow  loop  during  an 
experlMent  was  Monitored  by  periodically  Measuring  the  pressure 
drop  across  a  specified  p a*r  of  pressure  taps.  A  decrease  in  the 
Measured  pressure  drop  across  the~e  two  tars  would  indicate  that 
drag  reduction  due  to  additive  accunulatlon  was  occurring.  It 
should  be  noted  that  this  probleo  was  not  encountered  during 
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these  e iper iaents. 

The  additive  solution  use*  for  Mil  concentration  u«iur«- 
aents  m<  dyed  with  2.0  graas  per  liter  Flourescein  disodiua 
•alt.  Staples  aero  drown  froa  Mil  tape  (luring  the  injection  of 
the  dyed  solution  and  the  dye  corcentration  of  the  saaples  was 
aeasured.  The  saapling  rate  for  all  concentration  aeasureaonts 
was  SO  ml/min.  The  rationale  for  this  saapling  rate  it  discussed 
in  Append it  A. 
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RESULTS 


In  all  of  th#  experiments*  tho  Motor  temperature  in  tho 
channel  was  Mint  ai nod  at  24°C  and  tho  chonnol  flowroto  earn  840 
liter*  par  olnuta  yielding  a  Reynolds  number  of  17*800  based  on 
channel  width.  All  combination*  of  two  lovtla  for  each  of  tho 
independent  variable*  were  to«tod.  Th*  independent  variable* 
were  slot  angle*  slot  width*  injection  flowrate  and  injection 
concentration.  Tho  level*  of  these  variable*  and  the  combina¬ 
tion*  tested  are  susMsarited  in  Table  1. 


Table  1.  Experimental  conditions 

SLOT  ANGLE  SLOT  WIDTH  INJECTION  INJECTION 

FLOWRATE  CONCENTRATION 

< degroom)  <mm)  (al/ain)  <ppm) 


29  2.94  400  400 

-  "  -  100 

•  "  200  400 

"  -  -  too 

-  1.27  400  400 

"  -  -  100 

■  •  200  400 

■  '  "  100 

19  2.94  400  400 

•  -  -  too 

•  -  200  400 

•  •  -  100 

■  1.27  400  400 

•  "  100 

*  •  200  400 

•  "  *100 


Am  mtated  prev  ioumly,  drag  reduction  earn  deduced  free  prow 
mure  drop  aeaaureaentm  mode  in  the  channel.  Theme  aeaaureaentm 
provide  an  indirect  aeamurc  of  the  eall  mheer  mtreaa  in  the  chan¬ 
nel  and  therefore  a  measure  of  the  change  la  magnitude  of  th* 


12  - 


vlscout  drag  wNn  the  Additive  tolution  is  present.  For  e  fully 
developed  channel  flow*  the  average  well  sheer  stress  over  e 
given  length  of  the  chennel  is  proportionel  to  the  pressure  drop 
over  thet  length.  The  fully  developed  essueption  is  not  appllca- 
ble  in  the  vicinity  of  the  injection  slots  when  fluid  is  being 
injected*  nontheless  the  pressure  drops  still  yield  the  best 
estisets  of  the  viscous  dreg.  Hence*  in  this  study  dreg  reduc¬ 
tion  is  celculeted  fro* 


•Here  AP>4  is  the  pressure  drop  when  on  edditive  is  being 
injected  end  AP  is  the  pressure  drop  For  the  Fully  developed 
chennel  Floe  without  injection. 

Figures  3  through  6  shoe  the  variation  oF  drag  reduction  as 
a  Function  oF  d leans ion less  distance  <x*>  doenstrese  oF  the 
injection  slot.  The  distance  <x>  is  nor ea Used  eith  the  shear 
velocity  oF  the  channel  Floe  without  injection  and  the  It  In  seat  lc 
viscosity  oF  the  channel  eater .  The  level  oF  drag  reduction 
aoeeurod  between  two  taps  is  indicated  by  a  horizontal  line* 
spanning  the  distance  between  the  taps*  with  a  vertical  bar  at 
each  end.  The  vertical  line  e*vh  horizontal  bars  through  the 
point  indicates  the  uncertainty  For  a  tQX  conFidence  Interval. 
The  points  are  plotted  at  the  st re —else  location  side eg  between 
the  two  taps.  Each  Figure  cospares  the  porForsenco  *5  the  vari¬ 
ous  slot  geoeetries  For  a  given  injection  Flowrate  and  concentra¬ 
tion.  The  Four  Figures  cover  the  range  oF  Flowrates  and 


Comparison  of  various  slot  geometries  for  an  injection  flowrate  of  400  ml /min 
and  an  injection  concentration  equal  to  100  ppm. 
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Comparison  of  various  slot  geometries  for  an  Injection  flowrate  of  200  ml /ml 
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concentration*  t««t>d  in  this  study. 


The  dominant  characteristic*  of  Figures  3  through  4  it  the 
small  influence  of  the  slot  width  and  slot  angle  upon  the  result¬ 
ing  drag  reduction.  The  lack  of  geometric  dependence  shown  in 
Figures  3-4  could  be  the  result  of  a  constraint  Introduced  in  the 
choice  of  the  initial  levels  of  the  four  independent  variables. 
It  was  known  from  previous  experiments  that  when  She  injection 
momentum  f lux  through  a  0.  123  mm  slot  that  was  normal  to  the  flow 
direction  osceeded  approi imately  1/30  the  streamwise  momentum 
flux  through  the  viscous  sublager  that  the  Injected  fluid  would 
not  turn  and  flow  downstream  along  the  wall  but  would  jet  out 
into  the  main  flow.  In  the  interest  of  keeping  the  injected 
solutions  in  the  near-wall  region  for  as  long  at  possible*  the 
combinations  of  flowrates  and  slot  geometries  for  theso  initial 
eiperiments  were  chosen  to  that  the  normal  momentum  fluxes  were 
of  the  order  of  1/30  the  sublayer  streamwise  momentum  flux. 
PI  none  ion  low  normal  momentum  fluxes  <h+>  for  each  flowrate  - 
slot  geometry  combination  are  tabulated  in  Table  2.  8ince  slot 
width  it*  by  definition*  measured  in  the  streamwise  direction* 
normal  momentum  flux  is  dependent  on  slot  width  only  and  is 
independent  of  slot  angle.  Tho  values  of  momantum  flux  are  nor¬ 
malised  with  1/30  the  sublayer  momentum  flux. 
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Table  2.  Dimans ion lavs  normal  momentum  flux  of  injected  solutions. 


INJECTION  FLOWRATE 

SLOT  WIDTH 

H* 

(ml/min) 

(mm) 

400 

2.  94 

1.80 

400 

1.  27 

3.  60 

200 

2.  94 

0.  49 

200 

I.  27 

0.  90 

This  constraint  on  the 

injection  normal  momentum 

flux  did 

yield  injections  that  initially  remained  in  the  near-wall  region. 
This  was  confirmed  by  side  view  flow  visualisation  of  dyed  injec¬ 
tions.  For  all  of  the  geometries  studied*  the  additive  left  the 
injection  slot  as  a  thin  sheet  flowing  over  the  downstream  edge 
of  the  injector.  There  was  no  evidence  of  a  jet  catenating  from 
any  of  the  injectors  tested.  8ince  the  combinations  of  injection 
flowrates  and  slot  geometries  that  were  tested  did  not  cause  the 
additive  solutions  to  jet  away  from  the  wall*  the  only  influence 
of  slot  geometry  teas  to  modify  the  Initial  mixing  of  the  injected 
additive  with  the  water  flow. 

Figures  7  through  10  show  comparisons  of  the  perforaiance  of 
a  given  slot  when  tested  at  various  injection  flowrates  and  con¬ 
centrations.  The  four  figures  cover  the  range  of  slot  geometries 
tested.  It  can  be  seen  from  these  figures  that,  in  general*  the 
levels  of  drag  reduction  measured  Increase  with  mass  flow  rate  of 
polymer.  It  is  also  evident  that  drag  reduction  does  net 
increase  linearly  with  the  product  of  injection  flow  rate  end 
concentration. 

It  should  be  noted  that  Figure  10  shows  two  experiments 


Figure  8.  Comparison  of  various  injection  flowrates  and  concentrations  for  a  slot  angle 
of  15  degrees  and  a  slot  width  equal  to  2.54  run. 
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where  the  injection  flowrate  was  equal  to  400  ml/mln  for  an 
injection  concentration  of  400  pa.  This  coaparison  demonstrates 
the  rsproduceability  of  the  data. 

In  addition  to  the  drag  reduction  data*  aeasureaents  of 
additive  concentration  in  the  necr-wall  region  at  various  points 
downstreaa  of  the  injection  slot  were  aado  during  the  injection 
process  for  several  coabinatlons  of  the  Independent  variables. 
Figure  11  shows  a  typical  plot  of  near-wall  concentration*  as  a 
function  of  distance  downstreaa  of  the  injection  slots.  As 
before  the  distance  is  normalized  with  the  no-injectlon  channel 
shear  velocity  and  the  kinematic  viscosity  of  the  channel  water. 
Peak  values  of  concentration  for  various  esperlaental  conditions 
are  presented  in  Table  3*  a  more  comprehensive  compilation  of 
concentration  data  is  contained  in  Appendli  B. 

Of  Interest  here  is  the  magnitude  of  the  drag  reduction 
acMmvmd  For  the  levels  of  concentration  measured  from  • 
30*000  to  70*000.  Drag  reduction  levels  are  on  the  order  of  20X 
for  concentrations  less  than  1  ppm.  The  measured  concentrations 
(0.  A  ppm)  are  typical  of  the  fully  miied  concentration  for  this 
injection  flowrate  and  injection  concentration.  These  results 
deswnstrate  the  effectiveness  of  the  additive  once  it  has 
attained  a  drag  reducing  conf ormation  even  if  the  additive  con¬ 
centration  is  at  a  very  low  level. 
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Injection 

Flowrate 

(al/eln) 

Table  3.  Peak  concentrations 

Injection 

Concentration 
(P  P»> 

Concentration 

Slot 

angle* 23° 

Slot 

•idth-3.34 

OMB 

400 

400 

9.  1 

300 

400 

2.9 

Slot 

angl^lS0 

Slot 

•idth-3.34 

•• 

■ 

400 

400 

20.  0 

400 

100 

24.  33 

200 

400 

10.  9 

300 

100 

10.  3 

Slot 

angle-130 

Slot 

•ldth— 1 .27 

- 

400 

400 

14.  0 

400 

100 

4.  9 

300 

400 

9.  1 

300 

100 

3.3 
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REGRESSION  ANALYSIS 

Optimization  o#  the  additive  injection  process  was  conducted 
using  the  technique*  of  response  surface  methodology  based  upon 
factorial  eiporimentation.  To  implement  these  methods*  eaperi- 
ments  were  performed  over  a  small*  hut  significant*  rang*  of  tho 


Independent  variable*  and  the  response  (the  quantity  to  ho  SMilm- 
lied)  was  measured.  These  results  were  then  used  to  fit  a 


mathematical  model  of  tho  response  over  this  range  using  a  suit¬ 
able  regression  technique.  Once  the  mathematical  model  was 
determined*  the  importance  of  a  given  independent  variable  in 
determining  the  response  was  evaluated  statistically.  Also*  the 
direction  of  movement  in  the  space  of  the  independent  variables 
which  will  produce  the  largest  increase  in  the  rosponse  was 
determined.  The  latter  result  is  ef  prime  importance  in  the 
optimisation  process. 

Since  the  purpose  sf  this  study  is  te  melinite  tho  drag- 
reducing  perforsMnc*  of  the  additive  injection  process*  a  suit¬ 
able  measure  of  the  merit  of  a  given  esperimental  cenf iguratien 
was  devised.  This  measure, termed  the  merit  function*  was  defined 
as  the  area  under  the  drag  reduction  curve  from  the  point  where 
positive  drag  reduction  begins  <xq>  to  the  point  shore  the  injec¬ 
tion  process  ceases  to  he  effective.  This  definition  is 
equivalent  to  the  product  of  the  average  drag-reduction  and  the 
dimensionless  streemwis*  length  over  which  the  slat  injection  is 


effective. 
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With  th«  injection  of  on  additive  into  on  external  boundary 
layar*  one*  the  additive  diffuse*  fro*  the  near-wall  region  it  i* 
dispersed  into  the  free-«treo*  ond  lost  os  o  drog  reducer.  In  o 
chonnel  flow  the  additive  is  contained  ond  continually  mixes  bock 
into  the  near-wall  region  until  eventually*  there  is  only  o  very 
dilute  homogeneous  mixture  of  the  additive.  Once  this  homogene- 
ous  condition  is  reached*  any  drog  reduction  that  occurs  is  due 
only  to  tbe  presence  of  the  additive  ond  is  not  a  function  of  how 
the  additive  was  delivered  to  the  flow.  The  point  where  this 
homogeneous  region  hegln*  marks  where  the  diffusion  away  fro*  the 
near-wall  region  is  completed.  Since  only  the  effect  of  the 
injection  process  is  of  interest*  the  integration  for  the  merit 
function  stops  where  the  flow  becosM*  homogeneous.  Examination 
of  wall  concentration  and  drag  reduction  data  indicates  that  this 
point  is  in  the  vicinity  of  x*  •  90*000  for  the  condition*  oxom- 
ined  in  this  study.  In  practice  the  numerical  values  of  the 
merit  function  defined  above  are  quite  large  because  the  unit  for 
nen-diaensionelltlng  the  length  is  small.  Consequently  the  func¬ 
tion  <y>  used  in  this  report  was  multiplied  by  10  Thu*  the 
merit  function  Is 


-  30*000  . 

y  -  10  /  <XD*>dx 


Merit  function*  were  calculated  for  all  of  the  esperiSMnte 
and  the  results  wore  statistically  analysed  using  the  techniques 
of  response  surface  methodology.  As  discussed  by  Hunter  (1 *)« 
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the  first  step  in  the  analysis  is  to  define  a  sat  ef  naraaliiai 
coordinate  with  the  eriyin  located  in  the  cantar  ef  the  ran ye  af 
interest  i?f  the  Independent  varieties.  These  coordinates  aro 
defined  as: 

Xj  ■  (slot  angle  -  20°>/3° 

x2  -  (slot  width  -  1.88  m)/0.&3  mm 

x3  *  (injection  floerate  -  300  el/ein>/100  al/ein 

x4  m  (injection  concentration  -  230ppw>/130  ppa 

The  above  definitions  are  chosen  so  that  the  levels  of  the 
independent  variables  investiyated  becoee  plus  and  Minus  one  for 
all  of  the  independent  variables.  The  levels  and  ceabinations  of 
variables  tested  in  teras  of  these  new  variables  are  presented  in 
Table  4  alony  with  the  corresponding  oiperlaental  aerit  function 
<y>  for  each  coabinatlon. 

Since  only  two  levels  of  each  independent  variables  were 
tested*  the  aodel  chosen  is  linear  with  first  order  interactions. 

v  -  »B  »  V,  *  Va  ♦  *>3*3  *  w  biaxixa 

*  B13x1*3  *  bl4*lx4  *  b33xax3  *  b34*a*4  *  b 34 *3*4 

The  last  sic  teras  of  this  aodel  represent  interactions  between 
those  independent  variables.  Each  coefficient  represents  the 
relative  portion  of  the  total  variance  in  the  response  that  can 
be  attributed  to  a  ylvon  independent  variable.  A  variable  is  siy- 
nlflcant  only  if  the  variance  due  to  that  variable  is  yreater 
than  that  due  to  randoa  error.  It  should  be  noted  that  since  all 
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th*  *ip*rla*nt*  war*  net  r*plicat*d,  th*r*  it  no  way  to  separate 
randoe  *rror  f-oa  failure  of  th*  aodel  to  fit  th*  data.  There¬ 
fore  all  variance  away  froa  th*  fitted  aod*l  it  treated  at  randoa 
•rror. 


Table 

*1 

4.  Experlaental  and  aod*ll*d  value*  of  a*rit 

function  for  factorial  design 

*2  x3  *4  «  «• 

1 

1 

1 

820 

821 

827 

1 

1 

1 

-1 

329 

360 

499 

1 

1 

-1 

1 

714 

733 

674 

1 

1 

-1 

308 

339 

349 

1 

-1 

1 

1 

761 

766 

827 

1 

-1 

1 

**i 

498 

303 

499 

1 

-1 

-1 

1 

672 

680 

674 

1 

-1 

—1 

—1 

343 

284 

343 

-1 

1 

i 

743 

766 

827 

-1 

1 

1 

-1 

307 

303 

499 

-1 

1 

-1 

1 

713 

680 

674 

-1 

1 

—I 

—1 

226 

284 

343 

-1 

-1 

1 

1 

848 

821 

827 

*  -i 

-1 

-i 

399 

360 

499 

-1 

-1 

-1 

1 

730 

733 

674 

-1 

-1 

-1 

mm  f 

367 

339 

343 

Whan  this  aodel  was  fit  to  th#  experlaental  data  using  a 
least  squares  regression  and  the  statistical  significance  of  the 
est i sated  coefficients  was  tested,  the  only  coefficients  which, 
to  93X  confidence,  are  signif leant ly  different  froa  zero  were  b^, 
b4,  bja  and  b^.  The  coefficients  bja  and  are  significantly 
different  froa  zero  but  pi  eg  only  a  einor  part  in  predicting  the 
response;  the  bulk  of  the  variation  in  the  response  is  accounted 
for  by  the  coefficients  b3  and  b^.  The  final  fora  of  the  eodel 
is: 
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y^  -  586.13  +  76.75x3  ♦  164.25x4 
♦  27.75x^2  -  33.63X3X4 

Eu*in<tion  of  this  aodel  in  teras  of  the  itfinltiont  of  thi 
variable*  reveals  that  in  this  range  of  the  independent  vari¬ 
ables*  the  response  is  dominated  by  the  effects  of  injection 
flowrate  and  concentration.  The  fourth  and  fifth  toreis  in  the 
wodel  show  that  there  is  to  some  extent  an  Interaction  of  slot 
angle  and  slot  width*  and  injection  flowrate  and  concentration 
but  statistically  these  are  of  secondary  importance  to  the  Main 
effect*  due  to  concentration  and  flowrate. 

For  the  levels  of  the  independent  variables  tested  in  this 
study,  thw  rwsponsw  predicted  by  this  aodel  Cy^)  is  also 
presented  in  Table  4  along  with  predictions  calculated  using  just 
the  linear  portion  of  the  aodel. 

yj_  -  586.13  ♦  76. 75x3  ♦  164.25x4 

It  can  be  seen  that  the  prediction  of  the  five-tersi  aodel  is  in 
very  good  agreement  with  the  eiperiaental  results.  The  results 
predicted  by  the  three-term  linear  aodel  show  that  this  aodel 
accounts  for  the  bulk  of  the  variance  in  the  response.  This 
result  of  the  regression  analysis  confirsra  the  conclusions  based 
upon  inspection  of  the  drag  reduction  plots  in  Figures  3  through 
*. 

In  addition  to  identifying  the  relative  importance  of  the 
Independent  variables  to  the  merit  function*  the  regressien 


analysis  provide*  valuibli  information  for  tho  optimisation  pro** 
coot.  If  tho  osoct  mathematical  nature  of  the  rospono*  io  knoMn* 
ono  My  to  reach  on  optimum  rooponoo  Mould  bo  to  choose  o  start- 
inp  point  and  at  that  point  ealeulato  tho  fradiont  of  tha 
rooponoo  Mith  roopoct  to  tho  independent  variabloo.  This  diroe- 
tionai  dorivativo  give*  tho  direction  in  tho  doaain  of  tho 
indopondont  variable*  Mhlch  yields  tho  greatoot  incroaso  in  tho 
rooponoo.  Tho  rooponoo  is  thon  calculated  at  intarval*  along  tho 
path  prescribed  by  tho  gradient  until  a  poak  is  roachod.  Once 
this  poak  occurs*  tho  gradient  is  recalculated  and  tho  search 
takes  a  nee  diroction.  This  sequence  Mould  bo  ropoatod  until  a 
aaiiauo  is  roachod. 

This  procoduro  also  can  be  folloMod  in  an  esporlmontal 
optimization.  An  evaluation  of  tho  response  function  is  replaced 
by  an  esporlment  and  an  evaluation  of  tho  gradient  is  replaced  by 
the  sequence  of  factorial  esporlment*  described  previously.  Tho 
regression  analysis  is  the  kay  stop  in  evaluating  this  gradient. 
It  provides  a  mathematical  oiprossion  approslmating  the  local 
variation  of  tho  response. 

For  this  study*  tho  estimated  gradient  of  tho  response 
(merit  function)  Mhen  evaluated  at  tho  origin  of  tho  centered 
coordinates  yields  a  result  Mhich  is  indopondont  of  slot 
geometry.  Figure  IS  shoes  tho  diroction  of  steepest  ascent  in 
tho  ><3X4  plane  along  with  the  estimated  local  variation  in  the 
rooponoo  ee  predicted  by  the  linear  model  <y^>.  It  io  along  thio 
line  of  at oo poet  eocent  <k^*0.47x^>  that  future  expert mento 
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should  bo  conducted  in  order  to  approach  an  optimum  moot  effi¬ 
ciently. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  result*  of  this  studv  indicated  that  for  tho  rang*  of 
slot  geometries*  injection  flowrates  end  additive  concontrations 
examined*  the  performance*  depends  primarily  upon  tho  linear 
terms  of  injection  concentration  and  injection  flow  rate.  Tho 
product  of  injection  concentration  and  flow  rate  ao  wall  as  tho 
product  of  slot  angle  and  slot  width  woro  significant  but  loss 
iaportant.  The  statistical  trsataent  of  tho  results  also  yielded 
an  estlaate  of  the  local  gradient  which  in  turn  provides  tho 
direction  in  the  domain  of  the  independent  variables  to  aovo  in 
order  to  realise  the  largest  increaso  in  the  aorit  function.  The 
direction  indicated  involve*  increaso*  in  injection  flowrate  and 
additive  concentration  froa  tho  current  levels.  The  increase  in 
nondi mens  1  one  1  flowrate  <xg>  should  be  approximately  one  half  the 
increase  in  nondi mens ional  concentration  Cx^>. 

Since  the  path  for  future  experimentation  indicated  by  this 
study  will  involve  Increases  in  injection  flowrate#  hence  injec¬ 
tion  aoaentua  flux*  the  question  of  whether  the  injection  process 
is  inherently  independent  of  injector  geesHitry#  or  sms  rendered 
so  by  a  constraint  iaposed  on  noveal  aoaentua  flu*  from  the 
injector*  will  be  answered  during  the  course  of  further  experi¬ 
mentation. 

It  is  also  iaportant  to  net*  that  additive  concentrations 
less  than  1  ppa  yielded  drag  reduction  of  SOX.  It  is 
hypothesised  that  this  occurred  in  these  experiments  because  the 
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additives  had  attained  an  eicellent  conformation  for  drag  reduc¬ 
tion  by  the  time  they  reached  the  locations  where  this  result 


occurred. 
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Appendix  A  -  Concentration  Stapling  Technique 

Concentration  aeasureaents  were  aade  ae  part  of  this  study 
in  order  to  establish  a  data  baeo  for  future  investigation  of  the 
diffusion  characteristics  of  these  additives.  8ince  it  has  been 
shown  recently  in  our  laboratory  (1*2)  and  independently  verified 
(3)  that  these  additives  have  a  direct  effect  on  the  structure  of 
the  buffer  region  of  the  flow#  the  additive  concentration  in  this 
region  is  the  quantity  of  interest. 

All  concentration  aeasureaents  were  aade  using  saaples  of 
fluid  drawn  froa  the  near-wall  region  through  taps  located  flush 
with  the  surface  of  the  channel  wall.  The  taps  were  1.6  aa  in 
diaaeter  and  the  saaples  Here  collected  at  a  rate  of  20  al/aln 
over  a  tlae  period  of  one  to  one  and  one  half  ainutes.  This  sea- 
pi  ing  rate  is  the  lowest  practical  rate  which  can  be  aalntained 
using  the  present  facilities. 

Tba  rationale  for  using  this  s sapling  rate  <Qs>#  and  not  a 
higher  one  is  deaonstrated  in  Figure  A-l#  which  coapares  aeasured 
concentrations  for  two  saapllng  rates.  It  can  he  seen  froa  this 
plot  that  when  the  saapllng  rate  is  Increased  froa  20  al/ain  to 
40  al/ain  there  is  a  aarted  decrease  in  aeasured  concentration  in 
the  region  4,000  <  <  10000.  If  pooh  drag  reduction,  which 
occurs  in  this  region#  is  equivalent  to  a  peak  in  buffer  region 
concentration  then  this  decrease  in  aeasured  concentration  would 
indicate  that  the  fluid  swelled  at  40  al/ain  included  fluid  froa 
outside  the  buffer  region#  bonce  outside  the  area  of  interest. 


Variation  of  near-wail  concentration  (C*)  as  percent  of  injected  concentration 
ICj)  downstream  of  injection  slot  for  various  sampling  rates  (Qs). 
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Appendit  >  “  Concentration  Results 


«lot  angle: 

29  degrees 

•lot  width: 

2.  94  •• 

Injection  flowrate: 

200  «l/ein 

injection  concentration: 

400  ppe 

4*  —3 

x  xio  * 

concentration 

<PPe) 

2.  00 

11.  60 

10.  00 

0.  68 

22.  00 

0.  96 

94.00 

0.  36 

•lot  angle: 

29  degree* 

•lot  width: 

2.94m 

injection  flowrate: 

400  al/Mln 

injection  concentration:  400  ppM 

x  xio  * 

concentration 

(ppw) 

2.  00 

36.  47 

6.  00 

14.  90 

lO.  00 

9.  42 

14.  OO 

3.  32 

22.  00 

1.  13 

30.00 

.  9* 

38.  00 

.  9* 

46.  OO 

.  94 

94.  00 

.  44 

62.  00 

.  44 
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•lot  angle: 

•lot  width: 
injection  flowrate: 
injection  concentration: 


19  degree* 

1.  27  MR 
200  ml/min 
400  ppm 


x+xlO  ® 


c  one ent ret ion 
<P?«> 


2.  00 
10.  00 
22.  00 
94.  00 


20.  84 
1.  84 
.  64 
.  92 


•lot  angle: 

•lot  width: 
injection  flowrate 
injection  concentration: 


19  degree* 
1.  27  MR 
400  el/aln 
400  ppa 


x  xlO 


concentration 

(ppa) 


2.  00 
10.  00 
22.  00 
94.  OO 


96.  92 
9.  44 
1.  24 
0.  64 


•lot  angle: 
slot  width: 
injection  flowrate: 
injection  concentration: 


19  degree* 

1.  27  MR 
200  ml/min 

100  ppw 


X  XlO 


concentration 

(ppa> 


2.  00 
10.  00 
22.  OO 


3.  31 
.  23 
.  06 
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•lot  angle: 

•lot  width: 
injection  flowrate: 
injection  concentration: 


13  degrees 

1.  27  mm 
400  el /win 
100  ppm 


+  —3 

x  xlO  concentretion 

<PP»> 


2.  00 

3.  01 

10.  00 

.  38 

22.  00 

.  21 

34.  00 

.  11 

•lot  angle: 

13  degrees 

slot  width: 

2.  34  am 

injection  flowrate: 

200  al/aln 

injection  concentration: 

400  ppa 

4.  -3 

x  xlO  concentration 

(pp«> 

2.  00 

42.  48 

10.  00 

2.  28 

22.  00  t 

.  64 

34.  00 

.  48 

slot  angle: 

13  degrees 

slot  width: 

2.  34  sms 

injection  flowrate: 

400  al/ain 

injection  concentration: 

400  pps 

+  —3 ,  . 

x  xlO 

concentration 

(ppa> 

2.  00 

74.92 

10.  00 

3.  80 

22.  00 

1.  36 

34.  00 

.  72 

slot  angle: 

•lot  width: 
injection  flowrete: 
injection  concentration: 
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19  degrees 
2.  94  mm 

200  ml/min 
100  ppm 


concentration 

(ppa) 


2.  00  10.  29 

10.  00  .  69 

22.  00  .10 

94.00  .04 


•lot  angle:  19  degree* 

slot  width:  2.  94  mm 

injection  flowrate:  400  al/ain 

injection  concentration:  100  ppa 


+  -3 

x  xlO  concentration 

(ppa) 


2.  00 
10.  00 
22.  00 
94.  00 


24.  39 
1.  87 
1.  00 
.  26 
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Appendii  C  -  List  of  Publication*  and  Prosentaton 

The  following  list  sumsarixes  ttio  publications  and  preaenta- 
tions  during  tho  poriod  01  Harsh  1983  to  29  February  1984. 


1.  Luchik.  T.  8.  and  U.  0.  Tiadaraan.  Bursting  Ratos  in  Channol 
Flows  and  Drag-Reducing  Channel  Flowsi  Presented  at  the  Bya- 
posiua  on  Turbulence*  University  of  Hissouri-Rolla*  Sep¬ 
tember  26.  1983. 
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